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SOME  ASPECTS  OF  NON-EQUILIBRIUM  FLOWS 

SUMMARY 

In  this  paper  some  of  the  general  features  of  non-equilibri\am  flov 
are  discussed.  In  particular,  vibrational  relaxation  is  discussed  in 
detail.  This  case  is  somewhat  simpler  them  dissociation  and  ionization 
but  it  illustrates  some  of  the  main  new  featwes  of  non-equilibrium 
flow.  Those  aspects  of  two-dimensional  emd  axisymmetric  flow  behind 
shock  waves  are  examined  analytically  which  yield  significant  information 
without  requiring  numerical  solution  of  the  governing  eqxaations. 

The  thermodynamics  of  a  vibrational  relaxing  gas  is  discussed.  The 
conditions  for  simulating  flows  are  noted.  Crocco*  s  theorem  and  the 
characteristic  equations  are  derived.  Then  a  simple  method  of  obtaining 
the  initial  gradients  of  the  flow  variables  behind  a  shock  is  shown. 

These  gradients  axe  used  in  discussing  two  particular  flows.  An  exact 
solution  for  flow  over  a  cusped  body  is  obtained.  Flow  over  a  wedge 
near  the  tip  and  far  from  the  tip  is  considered.  It  is  foimd  that  far 
from  the  tip  a  boundary  layer  type  phenomenon  occurs. 
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pressvire  coefficient 

specific  heats  at  constant  pressiire,  volume  respectively 
2  2 

1  -  M  +  cot  X 
a 

operator  for  differentiation  along  characteristics 
internal  energy 

pi  (T^)  -  Ej  (T^  r  / 

coefficients  in  gradient  functions 

[®1  (fa'  -  \  ■"l' 

ds/Sp 

enthalpy 

cvirvature  of  shock  wave 
curvature  of  streamline 
displacement  of  shock.  Figure  7(c) 

Mach  number 

coordinate  along  normal  to  streamline 

pressure 

velocity 

gas  consteint  (per  moleciilar  weight  of  gas) 

radial  polar  coordinate 

entropy 

coordinate  along  streamline 

temperature 

time 

coordinate  in  free  stream  direction 
coordinate  normal  to  free  stream  direction 
Q  -  0 

shock  wave  angle 
ratio  of  specific  heats 

0  or  1  for  two-dimensional  or  axisyraraetric  flow  respectively 
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b 

Q 

o 

V 


vorticity 

coordinate  along  shock  wave 
direction  of  velocity  vector 
characteristic  vibrational  temperature 
p  -  9,  behind  shock  vave 
sin“-^  (1/mJ 


i 

JL 

P 

a 

T 

0 


coordinate  normal  to  shock  wave 
relaxation  distance  (equation  5*12) 
density 

distance  along  shock 
relaxation  time 
angular  polar  coordinate 
3  -  6 


a>  vorticity 

Subscripts: 

a  active  degrees  of  freedom  (translation  and  rotation) 

i  internal  degree  of  freedom  (vibration) 

e  equilibrium  conditions 

f  frozen  conditions 

t  stagnation  conditions 

oo  free  stream  conditions 
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1.  INTRODUCTION 


There  have  been  many  recent  publications  considering  the  flow  of 
a  gas  not  in  thermodynamic  equilibrlimi.  No  complete  bibliography  will 
be  given  here,  but  several  papers  pertinent  to  the  present  work  will  be 
mentioned. 


Probably  the  first  work  on  this  subject,  for  the  case  where  there 

are  shock  waves  in  the  flow,  is  that  of  Bethe  and  Teller^.  A  qualitative 

2 

discussion  of  flow  over  a  wedge  was  given  by  Ivey  and  Cline  in  which  the 

departure  from  thermodynamic  equilibriiom  was  due  to  vibrational  relaxation 

of  a  diatomic  gas.  (Some  quantitative  results  will  be  given  here  for 

this  case.)  An  approximation  to  the  hypersonic  flow  over  a  sphere,  where 

the  departure  from  equilibrium  is  due  to  dissociation,  was  considered  by 

5 

Freeman  .  Flows  with  small  disturbances  were  considered  by  Moore  and 

k 

Gibson  ,  and  in  particular,  the  flow  over  a  thin  wedge. 


In  a  pure  diatomic  gas  at  reasonable  tempera tiores,  three  processes 
can  cause  departure  from  equilibrium:  vibrational  excitation,  dissociation 
and  ionization.  If  the  relaxation  times  for  these  processes  are  suf¬ 
ficiently  different,  they  may  be  treated  separately.  The  magnitude  of  the 

■P'f*  "i"  n  -4“  Vk  -f*"!  /-VI  T  -?  r%  V.  4.1-^  _ _ 14.^  _ _ L 

\..j.  j.  uo  uiic  vcij.  j.ca.ua-co  lluc  uu  one  oiix  ee  ^X’uuetifciefc)  Xb  CJuixe  CLiiierenii 


The  gross  effects  are  directly  related  to  the  energy  necessary  to  excite 
the  new  degrees  of  freedom. 


In  this  paper  vibrational  relaxation  only  is  considered.  This  is 
somewhat  simpler  to  handle  than  dissociation  and  ionization,  but  should 
illustrate  some  of  the  main  features  of  non- equilibrium  flow.  Because  the 
smallest  non-equilibrium  effects  result  from  vibrational  excitation,  some 
of  the  results  obtained  here  differ  little  from  the  res\alts  of  calculations 
based  on  equilibrium  flow;  or  differences  are  significant  only  in  a 
negligibly  small  region  of  flow.  It  is  hoped  that  the  methods  used  here 
can  be  applied,  with  almost  equal  ease,  to  other  non-equilibrium  processes. 
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Only  certain  aspects  of  two-dimensional  and  axisyrametric  flows 
are  examined  in  this  paper;  these  will  not  involve  nvmerical  solution  of 
the  differential  equations.  After  a  discussion  of  the  appropriate  thermo¬ 
dynamics  of  a  vibrational  relaxing  diatomic  gas,  the  governing  equations 
are  transformed  to  yield  the  generalization  of  Crocco' s  theorem  which 
relates  the  entropy  change  to  the  vorticity.  Next  a  derivation  and  dis¬ 
cussion  of  the  characteristic  equations  are  given.  Using  the  appropriate 
shock  transition  relations  and  the  equations  of  motion,  the  streamwise 
gradients  of  the  flow  variables  are  obtained  employing  natural  coordinates 
(as  indicated  by  Sternberg'^^  for  equilibri\m  flows ^  These  prove  to  be 
quite  useful  in  discussing  two  particular  flows;  namely,  two  dimensional 
flow  over  (i)  a  cusped  body  which  supports  a  straight  shock  wave,  and 
(ii)  a  wedge* 

For  (i)  an  exact  solution  is  obtained  that  requires  only  two 
quadratures o  For  (ii)  the  flow  is  examined  near  the  tip  and  far  from  the 
tip*  At  the  tip  the  shock  curvature  and  the  wedge  pressure  gradient 
can  be  obtained  rather  simply  using  the  gradient  functions*  Far  from 
the  tip  the  flow  must  be  divided  into  two  regions*  To  a  first  approxi¬ 
mation  a  large  region  of  the  flow  is  the  equilibrium  wedge  flow,  but 
there  is  a  small  region  near  the  shock  wave  where  relaxation  is  im¬ 
portant-.  Mathematically  this  small  region  exhibits  a  boundary  layer  type 
phenomenon.  An  indication  is  given  of  a  means  of  obtaining  the  next 
approximat ion • 
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2,  THERMODYMMICS 


XJObtJ  ^Vw'Ut-L.-l-.y  UllV^  VJLA.  J.  A.  \wJi.  J>V^  ^  A.1X  .  Uixv^  VJ.V^tJV^A  X  — t-^LJX  XlXlll 


^V^u-X 


and  non-equilibrium  flows  results  from  differences  in  the  thermodynamic 
behavior  of  the  gas,  the  dynamic  aspects  being  the  same.  A  clear  explan¬ 
ation  of  an  appropriate  model  of  a  non-equilibrium  system  (and  the  necessary 

srsr  xr  ^ 

assumptions)  was  given  by  Wood  and  Kirkwood  .  The  same  model  will  be  used 
here. 

For  the  case  of  a  gas  subject  to  vibrational  relaxation  the  degrees 
of  freedom  are  divided  into  two  classes:  the  active  (translation  and 
rotation),  for  which  the  svibBcript  a  will  be  used,  and  the  internal  or 
inert  (vibration),  for  which  the  subscript  i  will  be  useda  It  is  assumed 
that  local  thermodynamic  equilibrium  exists  within  the  classes  but  not 
between  them.  The  rate  at  which  equilibriimi  is  approached  is  governed 
by  a  rate  equation  which  can  be  specified  in  terms  of  the  state  variables 
of  each  class.  There  are  other  assumptions,  to  be  mentioned  later,  which 
are  convenient  to  make.  The  model  presented^  is  intermediate  between  a 
macroscopic  sind  a  microscopic  description. 

From  the  assumption  that  the  Gibbs  relation  holds  for  the  “a*‘  class, 
it  ^TolloTf^s  th.a<t  th.e  entroj^y  ch.an.ge  of  th-e  "a"  class  is 

a 


where  the  perfect  gas  law  is  assumed 
p  =  p  R  T 


(2.1) 


It  is  assumed  that  the  1"  class  is  specified  by  its  temperat\ire, 
so  that  the  entropy  change  is 

n  /~i  .mn  /m  » 

^i  =  ^±^^1 

The  total  entropy  change  is  then 

AC*  _  .qo  1  ;jo 
UjC)  =  -T-  Ux:>, 

a  i 


=  c  d  T  /t  -  R  dp/p  +  dE./T. 
n  a  a  i  i 

“^a 


(2.2) 
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(2.5) 


The  rate  equation  assumed  is  the  linear  form 
dE^/dt  =  [e^(T^)  -  eJ  /t 

where  t  is  the  relaxation  time  and  (T^)  is  the  value  of  if 
equilibrium  existed  at  the  temperature  T^.  For  a  flow  the  derivative 
in  (2.5)  is  the  substantial  derivative. 

The  linear  form  of  the  rate  equation  presupposes  that  the  departure 
from  equilibrium  is  not  too  great.  It  is  exact  if  the  vibrators  are 
quantized  harmonic  oscillators.^  The  functional  relation  E^(T^),  for 
various  gases,  can  be  obtained  from  tables  of  the  properties  of  gases. ^ 
With  the  assumption  of  heirmonic  oscillators,  E. (T^)  can  be  calculated  from 

X  8* 

E^  (\)  =  R  ©y  [exp.  (ey\)  -  l|  ;  (2.*^) 

where  0  is  the  characteristic  vibrational  temperature.  This  form  will 
V 

be  used  when  some  specific  flows  are  calculated.  It  is  convenient, 
though  not  necessary,  to  carry  along  both  variables  E^  and  T^.  The 
functional  dependence  Ej^(T^),  for  the  case  of  harmonic  oscillators,  has 
the  same  form  as  E . (T  )  in  (2.4). 

X  Si 

For  future  reference  the  following  well  known  relations  are 
recorded.  For  equilibrium  flow  T^  =  T^  and 

c  -  c  =  R 


where 


Also 


c  =  c  +  c. 
up  i 
■a 


c  =  c  +  c^ 
V  Va  i 


c,  =  dE,/dT^. 


c  -  c  =  Ri 
p  V 

a 
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5.  EQUATIONS  OF  MOTION 


In  this  section  the  equations  of  motion  of  a  p’ore  diatomic 
gas  subject  to  vibrational  relaxation  will  be  discussed.  The  equation 
of  state  (2,1)  is  assumed.  By  a  minor  modification  the  resiilts  can 


be  extended  to  t*h6  csso  of  8^  niixtioro  of  diatomic  gases  vh6i*6  only  one 
constituent  relaxes.  It  is  necessary  to  know  the  effect  of  the  other 
constituents  on  this  one.  This  case  would  arise,  in  practice,  if  the 
relaxation  time  of  one  constituent  is  much  shorter  than  that  of  the 
others.  This  would  be  the  case  for  air  since  the  relaxation  time  of 
oxygen  is  about  one-fifteenth  that  of  nitrogen.  This  extension,  how¬ 
ever,  will  not  be  included  here. 


The  equations  will  be  written  in  natural  coordinates  (s,  n), 
where  s  and  n  are  distances  along  the  streamlines  and  their  orthogonal 
trajectories,  respectively.  The  form  of  the  equations  of  conservation 
of  momentum,  mass,  and  energy  do  not  change  because  of  relaxation;  they 
are 


p  q.  dq/ds  =  -  ^/ds 

(5.1) 

P  q  So/Ss  =  -  up/un 

(5.2) 

1  ho  1  da  dQ  e  sin  0 

(5.5) 

h  +  (qV2)  =  h^ 

(5.4) 

where  q  is  velocity,  9  is  the  velocity  direction,  e  =  0  or  1  for  two- 
dimensional  or  axlsvmmetrlc  flow,  resuectlvelv:  v  is  the  eno-rdlunte 

V  -  ^  - - - 1/  V _  —  •  V  W 

psrp8ridiGuls.r  to  tho  BucIs  of  syinmctryj  h  =  c  T  +  E,  +  p/p  is  the 
enthalpy;  and  h  Is  the  stagnation  enthalpy  (constant  along  streamlines), 
These  four  equations,  together  with  the  equation  of  state  (2,1)  Mid  the 
rate  equation,  give  a  complete  description  of  the  flow, 
coordinates  the  rate  equation  (2.3)  is 


n  /ha  — 


1?  /m  ^  =  1? 


A 


/  ‘Z  !=•  \ 
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The  relaxation  time,  t,  Is  a  function  of  pressure  and  temperature. 
In  the  foUoving  it  is  assumed  to  be  constant.  This  would  be  a  poor 
assumption  if  there  are  large  changes  in  pressure  and  temperature;  however, 
in  the  regions  of  flow  considered  changes  of  appreciable  magnitude  in 
these  quantities  do  not  occur. 


For  the  above  six  equations,  the  conditions  for  similar  flows  over 

geometrically  similar  bodies  are  easily  deduced.  In  addition  to  the 

requirements  of  fixed  Mach  number  and  ratio  of  specific  heats,  t  u/-^  must 

also  be  fixed  for  similar  flows.  Here  U  and  are  representative  velocity 

and  length,  respectively.  This  type  of  dimensionless  ratio  was  used  by 

5 

Freeman  .  Experimentally  the  appearance  of  the  new  length  scale  tU  meikes 
simxilation  more  difflcvilt.  Analytically  it  means  that  the  simple  "similar¬ 
ity  solutions”  (wedge  flow  and  Prandtl-Meyer  expansion)  no  longer  exist. 


g  _ 

A  useful  relation  for  equilibrium  flows  is  Urocco's  theorem.'  it 
relates  the  entropy  gradient  to  the  vorticity.  For  a  relaxing  gas  this 
relation  can  be  derived  as  follows.  The  expression  for  the  entropy 
change  (2.2)  can  also  be  written 


T^  dS  =  dh  -  dp/p  - 


d  E 


i 


Using  (5.^)  to  eliminate  h, 

T  dS  =  dh  -  q  dq  -  dp/  p  - 


d  E 


1 


From  this  the  directional  derivatives  of  S  in  the  s  and  n  directions 
are  evaluated.  After  elimination  of  the  pressure  by  use  of  the  momentun 
eqxiations  (3.1)  and  (3.2),  these  derivatives  are  expressible  as 

T^  as/ds  =  -  [1  -  (T^/Tj^^  bE^/bs  (3.6) 

T^  ^/^n  =  qC  -  &  SE^/^n  +  dh^/dn  (5.7) 

where  ^  is  the  vorticity, 

^  =  q  bQ/be  -  bq/bn 
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Note  that  dh^/ds  =0,  but  3h^/dn  =  0  only  for  isoenergetic  flows.  If 

the  flow  is  in  equilibritun,  T  -  H ,  (5.6)  gives  the  result  of  constant 

8*  1 

entropy  along  streamlines,  and  (3.7)  reduces  to  the  usual  form  of  Crocco’s 
theorem.  For  three-dimensional  \insteady  flow  the  vector  form  of  Crocco’s 
theorem  is 


S4/^t  -  4 


t 

X  0)  =  T  7  S  - 
a 


f 


f 

7E^ 


where  a>  is  the  vorticity  vector  and  vis  the  gradient  operator.  An 
equivalent  resiilt,  but  in  different  thermodynamic  variables,  was  given 

X.  10 

by  croer. 


It  is  convenient  to  eliminate  p  and  p  from  the  equations  of  motion 
and  introduce  S.  The  resulting  equations  for  isoenergetic  flow  are 

Se, 


,  P  .  r^n  r\Ci 


C  n  o  n  -Tk  O 

■Wi  ^ 


+  ( 


c  T 
Po  a 


) 


(5.8) 


P  . 


q  dO/ds  -  q  dq/^  =  T  ds/dn  + 


1  -  (T^T,)!  dE,/dn  (3.9) 


ds/ds  +  jl  -  (T^/T^^  dE^/ds  =  0 

/  r 

q  dEj^/ds  =  |1^(t^)  -  E^  /T 
%  '^a  \ 


(3.10) 

(5.11) 

(3.:i2) 


where 


c„  /(c_  R  T_), 


Thus  only  the  "frozen”  Mach  number  enters.  Note  that  (3.9)  and  (3.10)  are 

simply  (5.6)  and  (5.7)  while  (3.8)  comes  from  (2.1),  (5.1),  (3.5),  and 

the  expression  for  ds/ds.  The  unlcnowns  are  q,  9,  T  ,  E  ,  and  S. 

a  i 

It  is  easily  verified  by  the  standard  technique  that  the  above  partial 

differential  equations  are  hyperbolic  if  M  >1.  Solving  the  determinant 

a 

equation  for  the  four  characteristic  directions,  one  obtains  the  streamline 


15 


direction  coimted  twice  and  the  Mach  lines  based  on  M  ;  that  is,  the 

Mach  angle  is  ^  =  i\i\i  “■  (l/M  )•  The  equations  of  motion  in  eheiracteristic 

a 

form  axe  (5.10)^  and  the  follovlng  two  equations:  p 

4^4  sin  |i  E^(T  )  -  E. 

4.  TV  ~  TV  A  6  <1  sin  Q  sin  u  .  ^  Li'  a'  il 

cot  n  D .  q  +  q  D .  e  =  - * - >-  +  — 

+  +  y 


T  c  T 


-  cot  \i< 


S  +  [l  -  (T^/T^)]  eJ  /q 


(3.15) 


where  the  operators  B.  are  defined  by 

sec  |i  D .  =  d/Ss  +  tan  p.  d/dnj 


that  is,  indicates  differentiation  along  the  Mach  lines  inclined  at 
+  u  to  the^streamline •  The  appearance  of  the  "frozen"  Mach  lines  has 
been  noted  by  other  authors^ A  numerical  calculation^  using  the 
characteristic  form,  for  flow  over  a  wedge  is  in  progress  and  will  be 
reported  in  another  paper. 
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4.  GRADIEMTS  BEHIND  A  SHOCK  WAVE 


A  general  method  for  computing  the  flow  variable  gradients  behind 
a  shock  wave  in  a  two-dimensional  flow  was  developed  by  Thomas^^  for  a 
gas  not  subject  to  relaxation  effects.  A  simpler  method,  using  natural 

5 

coordinates,  was  indicated  by  Sternberg.  For  two-dimensional,  flow  the 
gradients  are  particilLarly  useful j  e.g.,  the  slope  of  the  streamlines  at 
shock  polar  points  in  the  hodograph  plane  can  be  obtained  from  them  (the 
Bo-cadled  "hedge-hog”  Introduced  by  Busemann).  It  is  found  that  the 
gradients  are  proportional  to  the  shock  wave  cvirvature. 


For  axisymmetric  flow  the  same  gradients  can  be  computed,  but  they 
are  linear  combinations  of  the  two  curvatures;  K^  in  the  merldlanal  and 
l/y  in  the  azimuthal  planes.  This  fact  limits  the  usefulness  of  the 

12 

gradients  for  this  case.  Extensive  tab\alations  have  been  made  recently 
of  the  coefficients  of  K^  and  l/y;  with  these  the  gradients  can  be  computed 
easily.  (The  two-dimensional  res-ult  is  obtained  by  setting  the  coefficient 
of  l/y  equal  to  zero.) 


For  a  uniform  flow  which  is  in  equilibrium  in  front  of  a  shock 
wave  but  not  behind  it,  the  gradients  will  depend  on  the  relaxation  time. 
Behind  the  shock  the  flow  variables  are  computed  from  the  shock  relations 

'PrfcT'  'P't ViAO'l~a«  *1  ca  v.^1  ci+ -I ci  ^  ^ 

1  ^ 
change  across  a  shock.  To  obtain  expressions  for  the  gradients,  the 

equations  of  motion  are  combined  with  the  shock  relations;  the  normal 

derivatives  are  expressed  in  terms  of  derivatives  along  the  shock  wave 

eind  along  the  streamline.  Referring  to  Figure  1, 

h/ba  =  cos  \  h/ha  +  sin  X 

where  o  is  distance  along  the  shock  wave,  p  is  the  shock  angle,  and 
X  -  p  -  0.  Since 

{h/hp)  = 

the  normal  derivative  can  be  written 


p/dn  =  ^  h/be  /  sin 
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The  normal  derivatives  in  (3.8)  and  (3.9)  are  eliminated  by  use  of  (4.1); 


=  0,  and  dq/dp,  dQ/dp,  etc.  can  be  computed  from  the  shock 
relations.  The  following  two  equations  then  resiolt: 

(M^‘'  -  l)  dq/qds  +  cot  X  dG/ds  =  (e/y)  sin  0  +  f ^  +  (k;^  k^/sin  X) 

\  '  •— / 


where 


cot  X  dq/ds  +  q  dG/ds  =  (fg  +  g  k^)/(q  sin  X)  (4.3) 


kj^  =  dq/dp,  kg  =  do/ dp 


=  (aEi/Ss)/=j,  - 

9. 

-P  ^  m 


/  I  c  r  T  q/ 
'  ^  p  a 


and  all  quantities  are  evaluated  immediately  behind  the  shock  waveo 
Solving  (4.2)  and  (4.3)  f’or  dO/Ss  and  dq/ds,  one  obtains 


dO/ds  =  +  (e/y)  Fg  +  (f^^  cot  x)/d 

(1/q.Qo)  dq/ds  =  F  +  (e/y)  F|^  -  q  D) 


(4.5) 


where 


rr  -1  *#  ^ \  / 1. 


q.  )  Q.  ^2 


^  ^  /  (m.  D 


a-?  1  ^ 

i-»  JL.XJ,  /V  y 


F^  =  (sin  0  cot  X)/d 


n 


F^  =  +  ^2  cot  X  -  q  kg  /(q^  D  sin  X) 


U 


Fj^  =  -  (q  sin  G)  /(q^  D) 


D  =  1  -  M  ^  +  cot^  X 

a 
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y  = 


The  F.  are  functions  of  M  ,  P,  and  y ,  and  have  been  tabulated  for 

1.4,  1.1^  —  Jt/2. 


00 


l6 


By  use  of  (5.1)#  (3.5)#  and  (5.4)  the  gradients  of  p#  p#  and 
can  be  obtained.  Since  f,  and  D  axe  positive,  the  effect  of  relaxation 

X 

is  to  increase  the  streamline  curvatvire  bQ/ha  and  decrease  the  velocity 
gradient  for  given  K  and  y.  For  reasonably  strong  shocks  at  ordinary 

V 

T__ ,  (T  )  can  be  neglected  compared  to  (T  ) . 
uu  j.  GO  X  a 


Even  for  two  dimensional  flows  it  is  seen  thdt  the  gradients  will 


not  be  simply  proportional  to  K  :  this  will  limit  their  usefulness.  Some 

-  -  -  y 

mrsTTkl  ■? /-vv\  T  T  Vizi  rrv»Q -I  a  01^+ 

W  J-.*.  J.  V/kJ  S/XXV^  XXV^A  V  WXXXA^  k#  • 


l?X 


L  , 


;  i  i* 


;; 


STiui’v) 

-  /.i  -1  r\  « 

IRL.  a::,  to. 


5.  AN  EXACT  SOLOTION 


As  mentioned  previously,  the  solution  for  flow  over  a  wedge  is 
no  longer  simple  when  the  gas  behind  the  shock  is  relaxing;  the  flow 
behind  the  shock  is  not  uniform  and  the  shock  is  not  straight.  This 
flow  will  be  considered  in  the  next  section.  Here  it  will  be  shown  that 
an  exact  solution  can  be  obtained  for  two-dimensional  flow  over  a  cusped 
body  (the  CTurvature  becoming  zero  asymptotically)  which  supports  a 
straight  shock  wave. 

Prom  (4.4),  with  e  =  0,  it  can  be  seen  that  if  the  shock  is  locally 

straight,  K  =  0,  the  curvature  of  the  streamline  at  the  shock  can  be 

w 

easily  calc\ilated.  Letting  K_  =  ^O/Ss,  a  convenient  non-dimensional 

t5 

measure  of  streamline  curvature  is  the  following  quantity: 


K  T  q  /e.  =  (R  q  cot  X)/c  q  D 

s  ^00 '  1  '  CD  ''  p  ^ 


(5.1) 


where 


a  =  (t^)  -  /ra 


(5.2) 


A  plot  of  the  expression  in  (5.1)  is  shown  in  Figure  2  for  a  rajige  of 

values  of  M  ajid  6.  For  the  case  of  pure  with  T  =  300°K,  p  =  60^ , 
CO  2  oo  /  f 

M  =6,  where  (2.4)  was  used  with  9  =  3556°K,  it  is  fovind  that  9  =  4l.l 
oo  '  '  V 

and  K  =5.5  fu  .  The  value  of  t  was  taken  from  Reference  8. 
s 

Since  (5.1)  has  no  physical  length  scale  in  it,  the  shock  wave  co\xld 

be  straight  over  any  finite  distance  and  K  would  be  the  same  everywhere. 

s 


.,0.  J 


Uonsiderin g 'Ciie  manner  m  wnicn  a  iTumenca-L  boiuuiun  ux  uhe  characteristic 

equations  would  proceed,  starting  from  a  straight  shock,  one  can  see  that 

ft1 1  the  flow  variables  depend  only  on  the  distance  along  the  normal  to 

the  shock.  A  "shock  oriented"  coordinate  system  6,  n  is  introduced,  as 

shown  in  Figure  5,  and  a  solution  with  flow  variables  independent  of  r\  is 

sought.  The  streamlines  will  be  parallel  cixrves  with  initial  curvatxire 

given  by  K  in  (5.1)  and  initial  slope  calculated  from  the  frozen  shock 
s 

relations.  If  the  streeunllnes  are  physically  meaningful,  any  one  of  them 
can  be  chosen  as  the  body.  Reflecting  this  curve  through  the  free  stream 
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direction  gives  a  pointed  cusped  body  which  supports  a  straight  shock 
wave.  Anticipating  a  result  from  the  next  chapter,  the  curvature  of 
the  body  should  approach  zero  and  the  angle  of  the  body  should  approach 
the  equilibrium  wedge  angle  appropriate  to  the  given  shock  angle.  These 
expectations  are  verified  by  the  solution  obtained  below. 

Eqmtions  (3*8)  to  (3.12)  are  transformed  to  the  t))  system  by 
S/^s  =  sin  (^  -  O)  3/^4  +  cos  (p  -  0)  h/hr\ 
b/bn  =  -  cos  (p  -  9)  b/bi  +  sin  (p  -  0) 

Setting  aJ.1  t)  derivatives  equal  to  zero  and  letting 

t  =  P  -  e. 


one  obtains  the  following  equations; 

(m/  -  1)  sin  ♦  II  +  q  =08  ♦  H  = 

^a 

q  sin  d9/d|  +  cos  ^|f  dq/d|  =  0 

as/^i  + 

q  sin  ijf 

=  T  +  +  (q2/2)  - 

•^a 

On  elimination  of  (5*^)  can  be  written 
dq/q  =  -  tan  \|f  d0 
Integrating, 

q  cos  ijA  =  q^  cos  (p  -  0^), 


(T^) 


/t 


(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


(5.8) 


where  the  subscript  f  denotes  values  computed  from  the  frozen  shock 
relations . 


Combining  (5.5)  with  the  differentiated  form  of  (5*7)  and 
eliminating  a  first  order  ordinary  differential  equation  is  obtained. 
The  result  of  integrating  this  is 

=  |T^  +  (q^^/R)  sin^  cot  (^  -  G^)  tan  t 


-  (q^^/R)  cos^  -  9^)  tan^  t 


(5.9) 


where  again  the  frozen  shock  relations  have  been  applied.  is  obtained 
from  (5.7),  with  the  aid  of  (5.8)  and  (5*9)* 

ili  =  hf  +  /R)tan^  \|f  -  (l/2)  sec^  ^  q^^  cos^  (p  -  0^) 

"  P  2  -I  (5.10) 

+  (q^  /r)  sin  (p  -  0^^  cot  (p  -  0^)  tan  t 


E, 


With  all  the  variables  expressed  in  terns  of  9,  (5.8)  can  now  be 
used  to  obtain  9  as  a  function  of  | .  Thus 

9 

'■  F(9)  d9 

0„ 


I  =  T  cos  (p  -  0^; 


F(0)  =  tan  t  (dE^/d0)/ 


(El  (T^)  - 


The  equation  of  the  body  is,  in  parametric  form  (parameter  Q-.  ), 

®b 


=  q^.  T  cos  (P  -  9  )  F(9)  d9 

D  J.  I 


A 


"'b  ^  "f  “  ®f^ 


F(9)/tan  >lrj  dO 
J 


The  above  integrations  are  carried  out  until  E . (T  )  is  equal 

X  a 

is,  equilibrium  is  attained.  In  principle,  and  become 


to  E^;  that 
infinite  for 
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this  condition;  actvially,  the  curvatvire  of  the  body  approaches  zero 
rather  quickly  and  the  gas  is,  for  all  practical,  purposes,  in  equilibrium 
at  finite  and 

A  representative  result  (obtained  by  numerical  Integration)  of  the 

above  formulas  is  shown  in  Figure  4  for  at  =  8,  p  =  65°,  =  500°K. 

Equations  (2.4)  was  used  to  obtain  E. (T  ),  with  ©  =  3338°K.  The  eingle 

of  the  body  changes  about  three  degrees  from  =  ^3.6°  to  ©^  =s  46.7  > 

the  latter  being  the  wedge  angle  that  corresponds  to  a  shock  angle  of 

65°  when  the  equilibrium’  shock  relations  are  used.  For  the  latter 

relations  see,  e.g..  Reference  15.  The  largest  change  in  angle  occurs 

when  ©^  and  M  are  close  to  the  conditions  for  shock  detachment  in 
f  00 

frozen  flow. 

Thus  ein  exact  calcxxlation  (with  two  quadratures)  of  the  flow  with 
a  straight  shock  wave  is  possible.  These  resvilts  can  be  used  to  check 
approximate  methods  of  computing  flows.  In  the  next  chapter  this  exact 
solution  will  be  used  to  investigate  a  portion  of  the  flow  far  from  the 
tip  of  a  wedge. 


It  is  desirable  to  define  some  measure  of  the  distance  over  which 

relaxation  effects  are  Important.  As  in  defining  the  thickness  of  a 

viscous  botmdary  layer,  there  is  some  arbitrariness  in  any  definition 

1). 

of  a  relaxation  distance.  Following  Moore  and  Gibson  ,  one  such 
definition  woixld  be 

■  CD 


^  KXJ 

I  f  - 


where  E^  is  the  final  equilibrium  value 
eq. 

by  E. (T  ),  the  value  of  the  integral  would 

1  81 

Integration  would  be  easier.  Calling  this 


of  E^.  If  E.  is  replaned 
eq. 

not  be  changed  much,  but  the 

relaxation  distance  6  , 

r 
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.00 

Jo  *- 


/m  /m  N 
/-iVXa/ 


d| 


=  IT  COB  (p  -  0^ 


)  ^  ^  I”  (dE^/d9)  tan  \|f  /  E^(T^) 


(5.12) 


vhere  0  1b  the  equilibrium  vedge  angle.  For  the  particular  example 
e 

cited  above 


=  .156  T  =  .500  q^ 


No  reason  was  given  for  the  particular  definition  of  a  relaxation  distaaice 
in  Reference  4;  it  is  euialogous  to  the  definition  of  displacement  thick¬ 
ness  of  a  viscous  houndBiry  layer®  In  the  next  chapter  a  different  method 
of  obtaining  a  distance  is  given  which  results  in  a  coefficient  of  t 
thirty  per  cent  higher  them  the  one  given  for 


6.  FLOW  OVjat  A  WJfiDGE 
2 

Several  years  ago  Ivey  and  Cline  discussed  the  two-dimensional 
flow  over  a  wedge.  This  was  mainly  a  qualitative  discussion.  However, 
it  was  pointed  out  that  at  the  tip  the  shock  wave  should  have  the  angle 
appropriate  to  the  frozen  shock  relations,  and  far  from  the  tip  the 
(smaller)  angle  appropriate  to  the  equilibrium  shock  relations*^.  Also 
the  shock  curve  is  concave  downward.  Here  some  further  quantitative 
results  for  this  flow  eire  shown. 


It  seems  most  natural  to  use  polar  coordinates  (r,  0)  with  the 
origin  at  the  tip  and  0  measured  from  the  free  stream  direction.  Trans- 
forming  Equations  (5.8)  to  (3.12),  with  e  ®  0,  to  polar  coordinates  gives 
2 

(Ma  -  l)(r  cos  a  hqj^r  +  sin  a  hq/M)  +  q(r  sin  (X  dQ/3r  -  cos  a  ho/hy) 
*  (q/Cp  Q!  SE^/^r  +  sin  a  3E^/dj^) 


q  (r  cos  a  do/^r  +  sin  a  So/^ji)  +  q(r  sin  a  'bqj'br  -  cos  a  3q/^) 


a  -  T  (r  sin  a  ^/3r  -  cos  a  ^/dji) 

M  *1  M  ^  I  f  "V*  C!  ^  VI  ^  ^  A.  / "N  \ 

^  -  uuB  u  oo/op) 

cos  CX  (j3fUJL  T  aJLn  u  OQ/Op)  = 

\l  -  (Ta/Tj^^  (r  cos  a  +  sin  a  SE^/S^) 

q(r  cos  a  ^E^/Sr  +  sin  a  =  (r/r)  Ej^(Ta)  - 


„  m  J.  1?  J,  frt^/o'S  —  V. 

''p/a  '  -  VH 


where  a  =  0  -  0 


(6,2) 


(fi-3) 

(6.4) 

(6-5) 


To  consider  the  flow  near  the  tip  a  dimensionless  variable,  r/xq  , 
is  Introduced.  Until  the  appropriate  boundary  conditions  for  wedge  flow 
are  Introduced,  the  following  discussion  applies  to  any  two-dimensional 
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flov  (e.g.,  the  modifications  to  Prandtl-Meyer  flow  over  a  corner). 
Assume  that  all  the  flow  variables  can  be  expanded  in  a  power  series 
in  r/rq.^;  then  it  can  be  shown  that,  to  within  0(r/Tq^^)  terms,  the 
flow  is  frozen.  The  form  of  these  series  would  be 


f(r,0)  =  f  (0)  +  (r/xa  )  f, + 

'  '  o''  '  '  OO  •  i''  • 


The  zeroth  order  terms  satisfy 


p 

-  l)  sin  ”  q._  cos  oc  hQ  /B0  =  0 

o 

cos  sin  =  0 

OS  /o0  =0 

/S0  =  0 

^o 


Equations  (6,6)  and  (6.7)  present  two  possibilities;  either 


2  2 

M  sin  a  =  1, 


(6.6) 

(6.7) 

(6.6) 

(6.9) 


which  is  appropriate  for  corner  flow,  or 

=  S©^/d0  =  0, 


Which,  after  applying  the  wedge  boundary  conditions,  gives  frozen  flow 
over  a  wedge.  To  within  an  error  0(r/Tq^),  the  shock  wave  is  straight. 


The  next  step  wotold  be  the  solution  of  the  differential  equations 
for  q^,  0,,  etc.  Althovigh  this  is  feasible,  it  will  not  be  done  because 
the  most  interesting  information  can  be  more  easily  obtained  from  the 
gradient  functions  of  Chapter  4.  Since,  for  the  wedge,  bQ/bs  =0,  the 
curvature  of  the  shock  at  the  tip  is,  from  (4,4), 


k;^  =  -  cot  X)/DP3^ 


The  frozen  values  of  the  variables  are  used  in  the  right  hand  side 
non-dimensional  form 


(6.10) 


K  tq  /e,  =  -  (Rq^  cot  X)/c^  q  DF, 
w  oo'  1  oo  P  1 

a 

A  plot  of  the  expression  in  (6,10)  is  shown  in  Figure  5  for  a  range  of 
M  and  S.  For  the  case  of  p\are  N„  with  T  =  300°K^  p  =  60°,  M  =6, 

Pqo  =  1  atm., 

0  =  41.1°,  K  =  4.1  ft"^ 
w 

where  (2.4)  was  used  with  0^  =  3556°  Kj  the  value  of  T  was  taken  from 

Reference  8,  and  F,  from  Reference  12.  For  the  same  conditions,  except 

M  =  10,  K  is  220  ft  ,  which  means  that  the  curvatirre  woiild  not  he 
oo  ^  w  ' 

detectable  in  any  ordinary  scale  of  experiment. 

With  K  known,  the  initial  velocity  gradient  along  the  wedge  can 
w 

be  computed  from  (4.5) »  and  then  the  pressure  gradient  from  (5<>l)»  A 
convenient  non-dimensional  form  for  the  pressure  coefficient  is 

(acyas)(Tq^/e^)  =  2(pq/p^q^)  [(R/°p^D)  -  (F5V‘^oo/®1^ 

where 

C  =  2(p  -  p  )/p  q  ^ 
p  •‘^oo^'^oo  00 

A  plot  of  the  expression  in  (6.11)  is  shown  in  Figure  6.  For  the  above 
example  (M^  =  6)  SC p/^s  =  l6,0  ft.  ^ 

Thus  the  shock  ciirvature  and  the  rate  of  change  of  the  flow 
variables  on  the  wedge  at  its  tip  can  be  obtained  rather  easily;  if  the 
details  of  the  flow  are  desired,  the  series  expansion  method  should  be 
used.  Unfortunately,  the  same  method,  i.e.  employing  the  gradient  fmctions, 
cannot  be  used  to  investigate  cone  flow.  The  use  of  natural  coordinates 
at  the  tip  of  a  cone  is  Impossible  because  of  the  nature  of  the  singularity 
there.  The  series  expansion  method  could  be  used,  but  the  differential 
equations  for  the  first  order  terms  would  be  much  more  difficult  to  soIvb 
them  those  of  the  wedge  case. 
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To  consider  the  flow  far  from  the  tip,  a  series  expansion  procedure 
is  again  adopted,  but  now  rq^/r  is  the  variable.  The  assumption  that 
the  flow  variables  are  analytic  in  this  variable  will  be  seen  to  be  not 
xmlformly  valid.  The  difference  between  this  case  and  the  previous  one 
can  be  seen  by  examining  (6.4).  When  ^/^q^  approaches  zero,  the  right 
hand  side  vanishes  eind  nothing  vinusual  happens.  When  xq^/r  approaches 
zero,  the  left  hand  side  vanishes  and  the  result  is  that  E ,  (T  )  =  E . ; 
that  is,  the  flow  is  in  equilibrium.  Since  derivatives  have  been  lost, 
a  bo\mdary  layer  type  phenomenon  (singxilar  perturbation)  must  be 
expected.  The  order  of  the  system  (6.1)  to  (6.4)  is  reduced  from  fo\ir 
to  three  by  this  limiting  process.  Also  the  resvilt  that  the  flow  is 
in  equllibriiun  is  inconsistent  with  the  boundary  condition  E^^  =  E^  at 
the  shock.  °° 

If  the  flow  variables  are  expanded  in  a  series  of  the  form 
f(r,  0)  =  f^°^  (0)  +  (x  <l^/r)  f^^^  (0)  +  .  .  .  . 

the  zeroth  order  terms  for  q  euid  6  will  satisfy  equations  of  the  same 
form  as  (6.6)  and  (6.7)  except  with  zero  superscripts  instead  of  sub¬ 
scripts  and  M  replaced  by  the  equilibrium  Mach  number.  Thus,  to 
within  an  error  0(xq^/r),  the  flow  is  equilibriiom  wedge  flow.  The 
shock  is  straight  and  the  flow  is  uniform  and  parallel  to  the  wedge. 

It  is  helpful  to  keep  in  mind  here  the  procediire  for  obtaining  the 
viscous  flow  over  a  body:  (i)  the  non -viscous  flow  Is  computed  ignoring 
one  boiindary  condition  (i.e.,  tangential  velocity  equals  zero),  (il)  the 
stmicture  of  the  apparent  discontinuity  (vortex  sheet)  at  the  body  is 
investigated  using  the  boxindary  layer  approximation,  (ill)  the  non- 
viscous  flow  over  the  body  plus  displacement  thickness  is  computed,  etc. 
For  the  present  problem  steps  analogous  to  (i)  and  (ii)  will  be  taken. 
Note  that  here  part  of  the  problem  is  to  determine  the  boundary;  that 
is,  the  shock  wave. 
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In  Figure  7(S’)  a  streamline  and  the  fictitious  "equilibrium  shock", 

predicted  by  the  first  term  of  the  series  expansion,  are  shovn  in  clashed 

lines.  The  transition  across  the  shock  is  governed  by  the  equilibrium 

shock  relations.  The  angle  of  the  equilibrium  shock  is  fixed  by  M  and 

oo 

the  vedge  angle.  Extended,  it  must  go  through  the  tip  in  order  to  satisfy 
the  conseirvation  of  mass.  The  actual  shock  (solid  line)  must  be  displaced 
upstream  of  the  equilibrium  shock,  but  its  curvatijre  must  vanish  to  within 
an  error  0(Tq^/r).  To  this  approximation  the  non-equilibrl\mi  flow  behind 
a  straight  shock  describes  the  flow  near  the  shock  but  far  from  the  tip 
(Figure  7(i>))*  This  flow  is  obtained  from  the  exact  solution  of  Chapter 
5.  Thus  the  "boundary  layer",  or  rapid  transition  region,  is  determinable. 
Compared  to  the  equllibrlvun  region  between  the  shcxjk  and  the  body  this 
transition  region  is  negligibly  small.  It  and  the  shock  merge  into  the 
equilibrium  shock. 

The  question  remains:  how  to  Join  the  two  regions;  or  rather,  how 
far  is  the  shock  displaced  from  the  equilibrium  shock?  This  is  answered 
by  again  invoking  the  conseirvation  of  mass.  A  transition  region  stream¬ 
line  must  asymptotically  approach  the  equilibrium  streamline  that 
originates  at  the  same  point  of  the  free  stream  flow  (see  Figure  7(b)), 

The  displacement  of  the  shock  from  the  equilibrium  shock  is  denoted  by  !>. 
From  the  condition  defining  it  is  easily  shown  that 

i  =  b  tan  j^an  -  tan  (p^  -  (6.12) 

where  6  is  the  equilibrium  shock  angle,  ©  is  the  wedge  angle,  and  b  is 
the  distance  (measured  perpendicular  to  the  shock)  between  the  shock  and 
the  asymptote  to  the  transition  streamline  (see  Figure  7(b)  or  5).  There 
is  no  convenient  expression  for  b,  but  it  is  easily  determined  graphically 
when  equations  (5.II)  axe  plotted.  For  the  conditions  of  the  example  of 
an  exact  solution  given  in  Chapter  5,  viz.,  p  =  65°,  ©  =  46.7*^,  M  =8, 
in  N2>  (6.12)  gives 

i  =  .177  Tq^  =  .0008  ft. 
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whereas  |  =  .136  q  t.  In  Figure  7(c)  the  two  regions  are  shown  according 

r  00 

to  this  approximation;  i.e.,  to  within  an  error  0(Tq^/r).  The  equilibrium 
shock  is  shown  extended  to  the  tip.  The  next  approximation,  with  ein  error 
0(Tq^/r)  ,  would  show  the  downstream  edge  of  the  transition  region  to  be 
curved. 

This  next  approximation  has  not  yet  been  worked  out.  The  following 
is  an  outline  of  how  this  might  be  done.  By  developing  gradient  functions 
for  equilibrium  flow,  the  curvature  of  the  equilibrium  shock  wave  could 
be  obtained.  Then  let  the  actual  shock  have  this  same  curvature.  For  this 
curved  shock,  a  shock  oriented  coordinate  system  could  be  used  to  attempt 
to  get  the  next  approximation  to  the  transition  region. 
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